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We study steady-state charge transfer across an interacting resonance-level model connected asym- 
metrically to two leads. For a linear energy dispersion relation of the leads, we calculate current- 
voltage characteristics of the model exactly employing the scattering Bethe-Ansatz of Mehta- Andrei 
and find symmetric transport showing the absence of diode effect. Next we study a lattice version 
of this model with a nonlinear dispersion for the leads using the Lippmann-Schwinger scattering 
theory. We find that the inclusion of nonlinearity in the leads' dispersion causes rectification for 
asymmetric junctions but does not rectify for asymmetric interactions and perfect junctions. The 
model in the latter case can be mapped into a model of a single noninteracting electron in higher 
dimensions. 

PACS numbers: 73.63. Kv, 72.10.Fk, 73.40.Ei 
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I. INTRODUCTION 



Rectification is considered as current asymmetry for 
the forward and the reverse bias. In the past years, rec- 
tification in nanoscale coherent systems has got a lot of 
interest. Current rectification by single asymmetric or- 
ganic molecules has been predicted [Ij and realized ex- 
perimentally 0, 11] . There are also several theoretical 
and experimental studies of charge 0, Q and spin ^ 
rectification in different mesoscopic semiconductor het- 
erostructures. The future application of modern molecu- 
lar electronics largely depends on the high-quality molec- 
ular rectifiers. Thus one needs to understand the basic 
mechanism of rectification by molecules, i.e., what are 
the necessary and sufficient conditions such that molec- 
ular junctions act as rectifier? The mechanism of rectifi- 
cation by molecules is highly debated [7l-[To|. 

The mechanism of current rectification in the original 
semiconductor p-n junction diodes or the Schottky diodes 
consisting of metal-semiconductor junction is a mismatch 
of band structures which creates a potential barrier that 
blocks the motion of carriers in one direction while al- 
lowing them to flow in the opposite direction. In fact, 
there are also microscopic studies in the recent past along 
this direction with different hybrid structures showing 
charge and thermal energy rectification. But one finds 
current asymmetry in molecular junctions or nanostruc- 
tures even for similar types of electrodes. Spatial asym- 
metry and nonlinear interaction between carriers are re- 
garded as the necessary conditions for charge rectification 
in these systems but it is still not clear what are the suf- 
ficient conditions for rectification [TTj . Quantum impu- 
rities are the simplest models for molecules or nanoscale 
heterostructures. Here we examine charge transport in 
a quantum impurity, namely, interacting resonance-level 
model (IRLM) connected asymmetrically to two leads. 
The equilibrium physics of the IRLM is well studied l^] , 
and recently the nonequilibrium transport in the IRLM 
has been received a lot of interest |13l - [l6j . We employ lin- 
ear and nonlinear energy dispersion of the leads. Surpris- 



ingly we find symmetric charge transport in the IRLM 
for linear dispersion of the leads even with different tun- 
neling junctions within the scattering Bethc-Ansatz ap- 
proach of Mehta and Andrei (MA) [13| . The inclusion of 
nonlinearity in the leads' dispersion causes rectification 
for asymmetric junctions but does not cause rectification 
for asymmetric interactions and perfectly transmitting 
junctions. Our model in the latter case can be viewed as 
a single noninteracting electron in the presence of elastic 
barriers in higher dimensions. 

Rectification in asymmetric impurities is essentially a 
nonlinear transport phenomenon. Study of charge trans- 
fer across out-of-equilibrium quantum imp urities has at- 
tracted much attention theoretically |13l - [20| as well as 
experimentally 21] for quite some time. The nonequilib- 
rium steady-state properties of quantum impurities can 
be investigated within the time-independent scattering 
approach. Recently MA [ll] have developed a nonper- 
turbative framework generalizing the equilibrium Bethe- 
Ansatz to compute steady-state properties of an IRLM 
connected symmetrically to left and right leads with a 
finite chemical-potential difference. They have employed 
a linear energy dispersion relation for the leads which 
is necessary for the application of the Bethe-Ansatz in 
their technique. Here first we apply the scattering Bethe 
ansatz framework to derive an exact expression for the 
charge current through the IRLM connected asymmetri- 
cally to the two leads with linear energy dispersion. For 
nonlinear dispersion of the leads, one can use different 
theoretical techniques such as the nonequilibrium Green's 
function formalism [17| or the ab initio first-principles 
calculations [To| . One needs to make many approxima- 
tions to apply these techniques for any interacting model, 
thus practically it is not possible to derive nonlinear 
transport in lattice models exactly. Recently we have 
studied nonequilibrium charge transport in quantum im- 
purities with sinusoidal dispersion of the leads using the 
Lippmann-Schwinger (LS) scattering theory [l^, [l^]. We 
here employ that method to study charge transfer in a 
lattice version of the IRLM with nonlinear dispersion of 
the leads and asymmetric junctions or interactions. Our 
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purpose in this paper is also to facilitate a critical dis- 
cussion on the above two methods investigating nonlinear 
transport in interacting quantum impurities. 



II. SCATTERING BETHE-ANSATZ FOR 
LINEAR DISPERSION 



The IRLM consists of a resonant level of energy ed con- 
nected to two leads via tunneling junctions of strengths 
ti and ^2 and Coulomb interaction U between the level 
and the leads. Then we apply standard manipulations 
for impurity models where we keep only the s angular 
modes around the impurity and linearize the bath spec- 
trum around the Fermi energies. The Hamiltonian of the 
system as chiral 1 — d field theories is given by 

H = ^ / dx ^pi{x)dipc.{x) + eddU+ ^{tiil;l{0)d 
+ t2i^l{0)d + H.c.) +U J2 iji{0)MO)d^d . (1) 



where we need to introduce a cut-off (bandwidth D) to 
make the model finite in the renormalized theory. Also 
we need to take same Fermi velocity for the both leads 
which we set to unity here. The current operator in 
this model is defined as / = j(ti'(/'i(0)d — t2ijl{0)d — 
H .c.) / {2\/2) . Once we compute the many-particle scat- 
tering eigenstate \'^)s for the asymmetric model, we can 
determine the steady-state current between the two leads 
by taking expectation of / in l^')^ 



(2) 



Now we transform the field operators 'tlji{x),^p2{x) of the 
two leads to a new set of even and odd field op erators 
4!e{x),ipo{x) using V'lla;) = (ti^leix) +t2'4)o{ x)) / ^Jtl +tl 
and ip2{x) = {t2'4'e{x) — ti'4'o{x)) / \/ 1\ + t\ . Under the 
transformation, the Hamiltonian in Eq.([T]) is decomposed 
into two parts of the even and the odd field operators. 
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He 
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He + Ho with (3) 
-I j dx ^Plix)dMx)+t{^pl{0)d + d^iJeiO)) 
+eddU+UiPliO)MO)d^d , 
-I I dx iblix)dM^) + UijliO)MO)d^d , 



with t = + ^i)/2. Then the Hamiltonian in Eq.® 
is exactly similar in the form of the Hamiltonian obtained 
by MA in their symmetric IRLM study after the trans- 
formation to a symmetric/antisymmetric basis. Thus the 
steps from here to evaluate the current are quite simi- 
lar to that of Ref . [l^l . First we calculate single-particle 
scattering states for the different boundary conditions 
(i.e., incoming electron from different leads) by solving 



the single-particle Schrodinger equation (here we have 
incorporated discontinuities at a; = following MA). We 
define the single-particle scattering states |l,p) for those 
with an incoming particle from lead 1. 



|1>-P> = I dx 
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+ {e^'^ - l)'j^e{x)4{x)) + "^^Six)d^] |0) 



where Sp — 2 a.Tctan[t^ /2{p — ed)] and Cp — t/{p — e^). 
We get the state \2,p) (those with an incoming particle 
from lead 2) from the state |1,J3) by interchanging si- 
multaneously the field operators ^\{x) and 'ip\{x) as well 
as the tunnelings and t2- The many- particle scatter- 
ing state is constructed from the single-particle scatter- 
ing eigenstates using the open Bethe-Ansatz framework 
[13| . For that we have to calculate two-particle S matrix 
by finding the two-particle scattering states for different 
boundary conditions. The linear dispersion in the leads 
gives the freedom to choose the two-particle S matrix be- 
tween all electrons to be the same; this helps to generalize 
the construction to many particle scattering state \^)s- 
Next one forms a Bethe-Ansatz basis of eigenstates for 
the noninteracting electrons in the leads. The nonequilib- 
rium boundary condition (namely, the different chemical 
potentials of the leads) has been incorporated in incom- 
ing particles' Bethe-Ansatz momenta {pj} in |\E')s which 
are determined by solving the Bethe-Ansatz equations. 

Now following the above stated prescriptions, we evalu- 
ate the many-particle scattering state \^)s and the corre- 
sponding Bethe-Ansatz momenta. Then we use Eq. ([2]) to 
find the steady-state current between the leads. Finally 
taking the thermodynamic limit at zero temperature, we 
get 
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dp[pi{p) - P2{p)] 



FiF 



2 



(p~ed)2 + (ri + r2)V4 



(4) 



i^/2, where the distribution functions 
1,2) satisfy following coupled equations 
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K{p, k) = 
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K{p,k)p0{k)dk{5) 
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TT {p + k-2ed)^ + ^{p-k)^ 



The Bethe momenta p in the lead a are filled from the 
lower cut-off {—D) up to fcg which is derived from. 
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Equation ffl correctly reproduces the result of Refs. 
|13| and [22| for the symmetric tunneling junctions, i.e., 
ti = ^2- One can solve the above coupled equations to 
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find the nonequihbrium distribution of the Bethe mo- 
menta using Wiener-Hopf method for J7 — >■ cx) or nu- 
merically for arbitrary U. As the Eq.([5|) determining 
the distribution functions of the Bethe momenta are in- 
dependent of the tunneling junctions ti and t2, we find 
from Eq.(j4|) that current between the leads remains the 
same if we interchange the tunneling junctions between 
the left and the right leads keeping the chemical poten- 
tial of the leads fixed. Thus the transport is symmetric 
for the forward and the reverse bias, i.e., for the chemical 
potentials [fi + V, fi] and -I- V], where V is the bias 
(we have set charge as unity everywhere). It shows that 
there is no diode effect or rectification in the current- 
voltage characteristics of this model even in the presence 
of the spatial asymmetry and the nonlinear interaction 
between electrons. 



III. LIPPMANN-SCHWINGER SCATTERING 
THEORY FOR NONLINEAR DISPERSION 

At this point the obvious question comes to our mind; 
what happens with the inclusion of nonlinearity in the 
dispersion relation of the leads? As we have discussed 
before we study the transport problem for nonlinear dis- 
persion using the LS scattering approach [l^, [2^. We 
here consider a lattice version of the IRLM; the Hamil- 
tonian is given by 

H ^ Ho + V , (6) 

OO 

where Hq = - '{clc^+i + 4+iCx) + ed"-o 

a; — — OO 

-(iic^iCo + ^24^1 + H.c.) , 
and V = Uin^iTiQ + U2noni , 

where nx — c\.Cx is the number operator at site x. 
denotes omission of a; = —1,0 from the summation. The 
Hamiltonian in Eq.® describes a resonant level of en- 
ergy erf at site being coupled with two noninteracting 
leads of spinless electrons modeled by one-dimensional 
tight-binding lattice. We set the lattice spacing and h 
to 1. Also we have taken an arbitrary strength for the 
tunneling junctions as well as the interaction between the 
resonant level and the left/right leads. 

The energy dispersion of a single particle with wave 
number k is given by Ek — — 2cosfc, where — tt < k < n. 
The wave function 4>k{x) for a particle incident on the 
resonant level from the left (with < fc < tt) or from the 
right (with — tt < A: < 0) can be found in terms of the 
tunneling ti, t2 and e^. The transmission probability |ifep 
turns out to be the same for wave numbers k and —k. We 
will find later that the two-particle current may not have 
this symmetry as a result of the interactions. We can de- 
termine the two-particle energy eigenstate for this model 
exactly [l9, 20]. The noninteracting two-particle energies 
and wave functions of Hq are given by i^k = E^^^ + Ef^^ 
and 0k(x) = (I>ki{xi)(l)k2ix2) - <l)kAx2)(f>k2ixi), where 



k = (fci, k2) and x = {xi,X2)- A scattering eigenstate (in 
the position basis) of the total Hamiltonian H is given 
by the LS equation 

^k(x) = 0k(x) + C/iif^Jx)V'k(-l,O) 

+ C/2ifBjx)^k(0, 1) (7) 

where if^Jx) = (x|G+(^k)| - 1,0) and i^^Jx) = 
(x|G+(£;k)|0, 1) with G+(£;k) = l/(Sk -Ho + le). The 
subscript k in the full scattering state V''k(x) in Eq.([7]) 
denotes the momenta of the incoming electrons. The 
momenta {k'l, fcj} of the scattered electrons can be quite 
different from the incident momenta {fci,fc2}, but must 
satisfy the total energy conservation after elastic scatter- 
ing which is given by cos/ci -I- cos/c2 — cos/cj -I- cosfc2. 
Now it is easy to prove the following properties from 
the above definitions, K'^ (— 1,0) = iC_Ek(0, 1) — Ko 
(say) and i^^jO,!) = XbJ-1,0) = (say). Then 
'0k(— 1,0) and V'k(0, 1) in Eq.([7|) are found in terms of 
these matrix elements. 

^""^ ' ' i-{u, + U2)Kq + umkI-kI) 

I (n. . K^U,M-i, 0) + (1 - UiKo)MO, 1) 
'^''^ ' ' 1 - (f/i + U2)Ko + C/1C/2 (K^ - Kl) 

Thus the two-particle scattering states are determined 
fully using the above ?/>k(— 1,0) and V'k(0, 1) in Eq.©. 
The many-particle scattering states for the nonlinear 
dispersion of the leads can be calculated within a two- 
particle scattering approximation. This is surely a per- 
turbative approach for many particles, but it can be jus- 
tified for weak interaction and/or weak tunneling with 
lower density of electrons in the leads. 

Now we calculate the steady-state current in this model 
at zero temperature. The current operator on the leads is 
defined by, jx = ~i{cl,Cx+i ~ c\^-^^Cx). First we find two- 
particle current j(fci,fc2) by taking expectation value of 
jx in the two-particle scattering state \4'-k) = |0k) + I'S'k) 
[from Eq.([7l)]. j(fci, ^2) = ji + jc + js, where current in 
the incident state is jj — {4>k\jx\4'\(.) — 2A/'(sinfci|ifeJ^ -|- 
sinfc2|iA:2p)j s-nd the contribution from the scattered 
wave functions are jc = (<^k| ja^lS'k) + (S'k|ja;|0k), and 
js = {S]i\jx\S]i). The normalization factor Af in ji will 
disappear in the many-particle current. We calculate the 
change in two-particle current, 6j{ki,k2) = jc + js nu- 
merically. We find that Sj(ki,k2) —Sj{—ki,—k2) if 
ti 7^ ^2 even for Ui ~ U2- This implies that the two- 
particle current change due to the interaction is quite 
different for the particles incident from the left as com- 
pare to the right even at the same energy. We have seen 
similar asymmetry in (25 j for another impurity model. 
This asymmetry in the two-particle current in the pres- 
ence of interactions and asymmetric junctions is the rea- 
son for rectification in the many-particle current, and the 
amount of rectification will be larger with increasing two- 
particle current asymmetry. So in Fig.(IT]), we plot the 
ratio Sj{—ki,—k2)/Sj{ki,k2) versus the energy of two 
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FIG. 1: (Color online) Plot of the ratio 

5j{~k\,~k2)/5j{ki,k2) versus the energy of two inci- 
dent electrons for ti = 0.2, t2 = 0.5 and ed = —0.5. 1^2] is 
kept fixed at 1.2 while |fci| is changed from 0.8 to 1.15 . 



incident electrons for the asymmetric junctions and dif- 
ferent values of [/i and U2- There is a large asymmetry 
in the two-particle current for a value of td correspond- 
ing to a two-particle resonance [l^l- Now we evaluate 
the many-particle current, j = ji + 6j. First we take 
the chemical potential of the left and the right lead re- 
spectively = — 2cos(fco + Sk) and ~ — 2cos(fco) 
with bias V = — l^R- Here fco + (— fco) is the 
highest occupied wave number of the left (right) lead. 
In the thermodynamic limit, the noninteracting current, 
ji = J'^^^'^'" ((ifc/27r)2sinfc|t/cp remains the same in mag- 
nitude when we reverse the bias. Thus we just need to 
find the many-particle current change 5j due to interac- 
tions to see rectification. Within the two-particle scat- 
tering approximation, 5j is given by 



'-^SmM. (8) 



(27r) 



We can write a similar expression for the Sj with the re- 
verse bias, i.e., fco (—[fco + Sk]) is the highest occupied 
wave number of the left (right) lead. It can be checked 
that the many-particle current change will be asymmet- 
ric, i.e., 6j{V) —Sj{—V), in the presence of the two- 
particle current asymmetry by inspecting the dj for the 
forward and the reverse bias. Thus we find here that 
an interacting quantum impurity acts as a rectifier for 
asymmetric coupling to the leads with nonlinear energy 
dispersion. 

Finally we simplify the lattice model in Eq.([6]) by con- 
sidering ti = t2 = I and ed = which corresponds to a 
perfectly transmitting impurity in the absence of inter- 
action. Now the single-particle state 0k (a;) = e''^^ and 
Kq, Ki are explicitly given by two-dimensional lattice 
Green's function g^^(x). Ko = 5^^(0,0) - 5+^(-l,l) 
and Ki = g+ (-1,-1) - g+ (-2,0) where gt{x) can 



be found in terms of the complete elliptic integrals |23{ . 
The two-particle current change in this simplified model 
has been derived analytically for arbitrary values of the 
Ui and C/2- Though the total current change Sj{ki, fc2) is 
same for the left /right lead but jc and js separately 
are different in the different leads. For fci,fc2 > 0, 
jc{x > 0) = 2T and jc{x < 0) = 0; for fci,fc2 < 0, 
jc{x > 0) = and jc{x < 0) = — 2T; and for fci,fc2 
opposite signs, jcix > 0) = T, and jc{x < 0) = — T, 
where 



T = 2Im[[/iV'k(-l,0)<^;;(-l,0) 
Similarly we find for fci , fc2 ^ 



U2MOA)cf>Uo,i)] .(9) 



js{x > 0) = -2[C/2|V.k(-l,0)p + UllMO, l)P]lm[ifo 



-U1U2 Im 
A*(2e-''?i 



dqi( 

■K 



A(2e 



^2iqi+iQ 



in sin Q 



-2iqi+iC, 



A-K sin Q 



(10) 



where A = '(/;k(-l, 0)V'k(0, 1) and - Eq^ = -2cosQ. 
While js{x < 0) is given by Eq.(fTO|) with Q being re- 
placed by —Q and the sign of the coefficient of Im[Xo] 
being positive. Interestingly now we find Sj{ki,k2) — 
—Sj{~ki,—k2) for arbitrary values of Ui and U2. Also 
6j{ki, k2) = for opposite signs of fci and fc2. So there is 
no rectification in the simplified version within the two- 
particle scattering approximation. The simplified model 
with two electrons can be mapped into a model of a non- 
interacting electron in two dimensions with two impurity 
sites of strength Ui and U2 at (—1,0) and (0,1). Sim- 
ilarly for three electrons the model can be viewed as a 
problem of a single electron in three dimensions with im- 
purity sites of strength Ui and U2 being placed on two 
infinite parallel lines, and this mapping can be extended 
for N electrons. It seems that there will not be any rec- 
tification in this simplified version even beyond the two- 
particle scattering approximation and it can be confirmed 
by considering three-particle scattering explicitly. 



IV. DISCUSSION 

To conclude wc find that the IRLM with different tun- 
nel junctions acts as a rectifier for nonlinear dispersion 
of the leads while it can not rectify for linear dispersion. 
The IRLM is integrable by the scattering Bethe-Ansatz 
for linear dispersion, but the Bethe-Ansatz technique is 
not applicable for nonlinear dispersion. The current in 
Eq. (|4]) is symmetric with respect to the bias as the distri- 
bution functions in Eq.([S]) are independent of the tunnel- 
ing junctions. Now the equations in Eq.([5]) are valid only 
for the energy of the resonance level (e^) being greater 
than the energy of the upper bounds (fcg and fcg) on the 
distribution in momenta in both the leads. The limi- 
tation on the validity of Eq.([5]) arises from the fact that 
the derivation of Eq. (O does not include the bound-state 
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contributions coming from the poles of the scattering ma- 
trix [22]. In fact, it can be shown that the two-particle 
scattering states contain a two-particle bound state and 
the corresponding two-particle current is asymmetric for 
asymmetric junctions even for a linear dispersion [2^ . In 
our study of lattice models using the LS scattering ap- 
proach fl^, [l^], we include the contributions from the 
two-particle bound states. Thus the symmetric trans- 
port in the IRLM with linear dispersion most probably 
arises only for the value of which does not include 
bound states. The scattering Bethe-Ansatz technique of 
Mehta-Andrei is incomplete as it does not capture the 
contributions from the bound states. 

The Bethe-Ansatz momenta {pj} in the scattering 
Bethe-Ansatz method are determined using periodic 
boundary conditions in an auxiliary algebraic Bethe- 
Ansatz problem and it has been claimed that for infi- 
nite periodicity these momenta will coincide with those 
of the physical systems. But it is not clear what is the 
mechanism of dissipation (or exchange of energy) of the 
scattered particles in the leads within such an approach. 
On the contrary, in the technique used here for the lat- 
tice model, it has been assumed as the original Landauer- 
Biittiker (LB) scattering approach that all the dissipation 
occurs in the reservoirs connected to the leads. One can 
show an one-to-one connection between the LS scatter- 
ing theory and the LB scattering approach. The present 
technique based on the LS scattering theory has been 
successfully applied to study transport in different inter- 
acting mesoscopic systems [19, 20, 25]. Also the many- 
particle current within this technique merges with the 
current derived from the LB approach for noninteracting 
systems [ll|. Though there are a few numerical studies 
such as Rcf. 16] using the density matrix renormaliza- 
tion group (DMRG), to our best knowledge this is the 
first analytical study of nonlinear transport in symmet- 
ric/asymmetric IRLM for nonlinear dispersion. 



Recently it has been shown in Ref.fl^ that thermal 
transport is symmetric/asymmetric with respect to tem- 
perature difference in mono-mode mediated energy ex- 
change between two metals for the linear/nonlinear dis- 
persion of the metals. But in that study the two metallic 
leads are coupled by a single harmonic modes, i.e., linear 
interaction. Here we emphasize that in the IRLM due to 
local nonlinear interaction particles can exchange energy 
after scattering. Thus we expect to find rectification in 
the IRLM with asymmetric junctions for both the linear 
and the nonlinear dispersion of the leads. Finally, we 
can understand physically the mechanism behind recti- 
fication in the lattice model with asymmetric tunneling. 
Let we consider scattering of two electrons from one lead 
to another mediated by nonlinear interactions at impu- 
rity site. Two electrons with incident momenta {ki,k2} 
can scatter into different channels with momenta {k'ljk!^} 
satisfying the total energy conservation. The asymmetry 
in tunnel junctions creates a difference in the redistribu- 
tion of momenta after scattering from the interactions for 
electrons coming from the left or the right leads. Thus 
we find the two-particle current asymmetry with asym- 
metric junctions. For many particles if a finite bias is 
applied across the impurity then the asymmetry in the 
redistribution of momenta persists and that generates the 
asymmetry in current. 
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